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Abstract 

Quantum Zakharov equations are obtained to describe the nonlinear interaction between quan- 
tum Langmuir waves and quantum ion- acoustic waves. These quantum Zakharov equations are 
applied to two model cases, namely the four-wave interaction and the decay instability. In the 
case of the four-wave instability, sufficiently large quantum effects tend to suppress the instability. 
For the decay instability, the quantum Zakharov equations lead to results similar to those of the 
classical decay instability except for quantum correction terms in the dispersion relations. Some 
considerations regarding the nonlinear aspects of the quantum Zakharov equations are also offered. 
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I. INTRODUCTION 



The importance of quantum effects in ultra-small electronic devices,— in dense astro- 
physical plasma systems 2 and in laser plasmas* have produced an increasing interest on 
the investigation of the quantum counterpart of some of the classical plasma physics phe- 
nomena. For instance, quantum plasma echoes,— the expansion of a quantum electron gas 
into vacuum,— the quantum two and three stream instabilities^ and the quantum Landau 
damping 7 have been the subject of recent investigations. Also, quantum methods like the 
Wigner-Moyal transform have been used in the treatment of the Landau damping of classical 
partially incoherent Langmuir waves.— Quantum models like the Schrodinger-Poisson system 
have also been used, through the correspondence principle, for the numerical simulation of 
the Vlasov-Poisson system** 

In this context, a mathematical formulation, based on the quantum hydrodynamical 
model for charged particle systems , 10 ! 11 was introduced to study the quantum version of 
the low frequency ion-acoustic waves.— In micro-electronics, the quantum hydrodynamical 
model describes** negative differential resistance associated to resonant tunnelling diodes. 
It can also model ultra-small high-electron-mobility transistors-^ The quantum hydrody- 
namical model for charged particle systems was also successfully used for the description of 
quantum dissipation,— under the same closure hypothesis as that adopted by Manfredi and 
Haas. 11 In the case of quantum ion- acoustic waves,— several features of pure quantum origin 
were observed for the linear, weakly nonlinear and fully nonlinear waves. The linear quan- 
tum ion-acoustic waves are described by a dispersion relation which tends to the classical 
dispersion relation as quantum effects goes to zero, in accordance with the correspondence 
principle. The weakly nonlinear quantum ion-acoustic waves are described by a modified 
Korteweg-de Vries equation depending on a scaled H parameter. Finally, the fully nonlin- 
ear quantum ion- acoustic waves can have a coherent, periodic pattern, not present in the 
classical case. This points to the intrinsically more reversible features of quantum plasmas, 
as seen for instance in quantum echoes* and coherent patterns in the quantum two stream 
instabilities*. 

The purpose of the present paper is to continue this investigation by studying the non- 
linear coupling between the quantum ion-acoustic waves and the quantum Langmuir waves. 
At the classical level, a set of coupled nonlinear wave equations describing the interaction 
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between high frequency Langmuir waves and low frequency ion-acoustic waves was first de- 
rived by Zakharov.— Since then, this system have been the subject of a large number of 
studies.— In one-dimension, the Zakharov equations can be written (in normalized units) as 



where E is the envelope of the high frequency electric field and n is the plasma density 
measured from its equilibrium value. The system (0121) can be derived from a hydrodynamic 
description of the plasm a 1 ? 1 ^ by distinguishing two different time scales, the slow time scale 
of the ions and the fast time scale of the electrons. The low mobility of the ions as compared 
to that of the electrons justifies this kind of treatment. Since the Landau damping of the 
Langmuir waves is neglected in the fluid description, the model (012)) is restricted by the 
condition k <C ko-, where k is the wavenumber and kp is the Debye wavenumber. Also, a 
weak turbulence condition is to be satisfied.-^ 

In this paper, modified Zakharov equations are obtained by use of a quantum fluid ap- 
proach. Specifically, we assume a two species, one-dimensional quantum plasma in the 
electrostatic approximation. Pressure effects are neglected for the ions whereas the elec- 
trons are described by an isothermal equation of state. Contrary to the quantum degenerate 
case,— the present model is more suitable to investigate the classical limit h — > O.We do not 
include quantum statistical effects in the present investigation, and therefore, only quantum 
diffraction effects, responsible e.g. for tunnelling, are taken into account. 

The paper is organized as follows. In Section II we write the quantum hydrodynamic 
model for a two-species plasma and derive the Langmuir mode for quantum plasmas. In 
Section III, we obtain the quantum Zakharov system through a procedure similar to the 
classical one where a two-time scale formalism is used. In Section IV we study the influence 
of quantum effects in two relevant parametric instabilities: the decay instability and the 
four-wave instability. Section VI is devoted to a preliminary discussion of the nonlinear 
aspects of the problem and some of the remaining open questions. Section V is reserved to 
the conclusions. 
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II. QUANTUM LANGMUIR WAVES 



Before considering the nonlinear coupling between ion-acoustic and Langmuir waves, we 
examine the linear stability analysis of the Langmuir waves in the quantum regime.— For 
this purpose we consider a one-dimensional quantum system, composed of electrons and 
singly charged ions. The quantum hydrodynamic equations in this case become^ 2 - 

dne d(n e u e ) = Q 
dt dx 
drij + djniUj) = Q 

dt dx 

du e ^ du e = e_ E 1 dP e | h 2 d ( d 2 ^fn~ e /dx 2 \ 

dt 6 dx m e m e n e dx 2m 2 dx \ ^frT e J ' 
duj ^ dm _ _e_ E 
dt 1 dx rrii 
dE e . , 

— = — (ni-n e ), (7) 
dx e 

where E is the electric field, P e is the electron pressure, and n e , rij, u e , Ui, m e and rrii represent 
the density, fluid velocity and mass of electrons (e) and ions (i) , respectively. In addition, e 
and h are the vacuum dielectric and the scaled Planck's constants. Since we are interested 
in high frequency waves, the ion density n« can be assumed constant, at this stage. The 
pressure P e is obtained from an equation of state for the electrons, which basically depends 
on the thermodynamic properties of the system. In the present investigation, we consider 
the isothermal equation of state P e = KBn e T e , where T e is the electrons' temperature and Kb 
is the Boltzmann's constant. In view of their large mass, ions are treated classically. Also, 
in a first approximation, we consider cold, zero temperature ions. The Bohm potential term 
proportional to h 2 in (jSJ) is responsible for negative differential resistance in semiconductor 
devices^ and is associated to tunnelling. 

Linearization of the electron equations ©, (jSJ) and around the homogeneous equilib- 
rium n e = rii = uq, u e = and E = produces the following dispersion relation: 

. 2 = u J 2 e+ v 2 e k 2 + ^. (8) 

In Eq. (jHJ), to is the wave frequency, k is the wavenumber, to e = (n e 2 /m e e ) 1 / 2 is the electron 
plasma frequency and v e = (K£T e /m e ) 1//2 is the electron thermal velocity. Notice that both 
classical and quantum modes can be obtained from Eq. (jHJ). In fact, the classical limit h — » 
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gives the classical Langmuir wave dispersion relation.— According to Eq. (|SJ), the frequency 
lu is always real, and instability (or damping) of this wave cannot be observed. The main 
purpose of this paper is to obtain a model describing the exchange of energy between the 
quantum Langmuir modes shown above and the recently found quantum ion-acoustic plasma 
modes.— 



III. QUANTUM ZAKHAROV EQUATIONS 

In order to obtain the set of equations describing the nonlinear interaction between Lang- 
muir waves and ion-acoustic waves, in the quantum regime, we follow the derivation orig- 
inally made by Zakharov.— A general discussion of the validity of the Zakharov equations 
can be found in the review paper by Thornhill and ter Haar.— 

We first separate all fluid variables into high frequency (subscript h) and low frequency 
(subscript I) components, 

n e (x,t) = n + ni(x,t) + n h (x,t) , (9) 
rii(x,t) = n + ni(x,t) , (10) 
u e (x,t) = ui(x,t) + u h (x,t) , (11) 
Ui(x,t) = u t (x,t) , (12) 
E(x,t) = Ei(x,t) + E h (x,t) . (13) 

Notice that the high frequency portions of the ion quantities [Eqs. ffTUj) and (JTJJ)] were 
ignored due to the large ion mass. Also, from the very beginning we assume that departures 
from the quasi-neutral regime (n, ~ n e and ~ u e ) are provided only by the high frequency 
components of the electrons motion. The high frequency term of the electric field can also 
be written as 

E h {x, t) = ~E{x, t)e-^ + c.c. , (14) 

where E(x,t) is the slowly varying envelope of the high frequency term and c.c. refer to 
complex conjugate. Using the high frequency components of Eqs. (J^HZJ), we obtain, by the 
same procedure used in the classical case, 17 

.dE lv 2 e d 2 E h 2 d i E_u e n l ~ 
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where the term \d 2 E\ <C \oj e d t E\ has been neglected. Equation (|T3j) describes the evolution 
of the slowly varying amplitude E, as defined in ()14|) . 

We next proceed with the derivation of the equation for the low frequency part, ni, of 
the departure from the equilibrium density uq. After averaging over the fast time scale, we 
get a set of equations describing the low frequency part of the electron continuity equation, 
electron force equation and ion force equation, 

^ + n °^ = ' ^ (16) 

dui e g K B T e drii h 2 d 3 rii e 2 d\E\ 2 
dt m e 1 nQm e dx Amino dx 3 Am 2 e uj 2 dx 

' El = 0. 18 

at mi 

Convective terms were disregarded in view of a weak Langmuir turbulence assumption, as 
detailed by Thornhill and ter Haar.— Eliminating ui and E[ from Eqs. (j!6H18|) and assuming 
m e /mi < 1, we obtain 

d 2 m c2 d 2 m | h 2 d A m _ e d 2 \E\ 2 
dt 2 s dx 2 Amim e dx A Ami dx 2 ' 

where c s = (KBT e /mi) 1 / 2 is the ion-acoustic velocity. We call Eqs. (|15j) and (JT§j) the quantum 
Zakharov equations. 

For the following analysis, it is most convenient to normalize Eqs. (II 5|) and f)19|) . Nor- 
malized quantities are expressed as 

i = 2J— r) t = 2— u e t, 20 
V m, X P mj 



n = -— , E=J — - — —E , 21 

Am e no y 16 m e no T e 

where A e is the electron Debye length. In addition to (j2()H21|) . we introduce the dimensionless 
quantum parameter 

H =—k> ( 22 ) 
where Ui = (n e 2 /m^eo) 1 / 2 is the ion plasma frequency. The resulting system reads (we 
dropped bars for the sake of simplicity) 

dE d 2 E TT2 d 4 E 

'm + w- H w =nE ' (23) 



d 2 n d 2 n „o<9 4 n d 2 \E 
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+ tf7£u=^br- (24) 



dt 2 dx 2 dx 4 dx 2 
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The quantum parameter H given in (|22|) expresses the ratio between the ion plasmon energy 
and the electron thermal energy. This is to be compared with the dimensionless parameter 
characterizing quantum effects in the two-stream quantum instability— given by the ratio 
between electron plasmon and thermal energies. Here, the presence of ion-acoustic modes 
forces the appearance of ionic (inertia) parameters. Notice that for dense plasmas, 2 with 
particle density about 10 25 — 10 32 m~ 3 and temperature about 10 5 — 10 7 K, the parameter H 
is not irrelevant at all. For a completely ionized hydrogen plasma in these ranges of densities 
and temperatures, H typically goes from negligible values of order 10 -5 up to values of order 
unity. The presence of large values of H points to the possible experimental manifestation of 
quantum effects in the coupling between Langmuir and ion-acoustic modes in dense plasmas, 
particularly in astrophysical plasmas. 

In the next section, the model (j23H24|) is used to investigate two parametric instabilities 
extensively studied in classical plasmas: the decay instability and the four-wave interaction. 



IV. PARAMETRIC INSTABILITIES 



A. Decay instability 

Following strictly the treatment for the classical decay instability, 17 consider the proposal 

E = E e iikox -" ot) + £ 1 (t)e i(fcl ™* ) , (25) 
n = n%{t) cos(Kx - Sit) , (26) 

now for the quantum Zakharov equations (j23H24j) . where Ei(t) and n\(t) are first-order 
quantities, and 

cu = k 2 + HX, (27) 

Ul = k\ + R 2 k\, (28) 

Q 2 = K 2 + H 2 K A . (29) 

Notice that Eq. (|29j) is identical to the quantum dispersion relation obtained by Haas et 
alJZ (Section V) except for the value of H which, due to the use of a quantum equation of 
state, has a different definition there. 

Also, there are the usual matching conditions 

k = ki + K , ojq = Ui + Q , (30) 
7 



corresponding to momentum and energy conservation respectively. These conditions de- 
scribe the decay of one quantum Langmuir wave, with dispersion relation (J27)l . into other 
quantum Langmuir wave, with dispersion relation (|28|). and a quantum ion-acoustic wave, 
with dispersion relation (J29j) . 

Linearizing the quantum Zakharov equations, we obtain 

iEid = ^7p-(e+ + e_)e , (31) 
- zfirii + K 2 E Q E{) e + + (y + iQri i + K 2 E*E^j e_ = , (32) 

where use has been made of the notation 

eo,i = exp[i(k ^x — u>o^t)] , e± = exp[±i(Kx — Qt)] . (33) 

The resonant part [satisfying (jSOJ)] of (|3*T|) gives 



2? • 

ni = —E x , (34) 



while (J3~2*j) gives 



n 



- + zQrix + i^E*^ = . (35) 
For £i = iuEi, elimination of ni leaves us with 

cu 3 + 2Qu; 2 + JT 2 |^ | 2 = 0, (36) 

which is formally identical to the dispersion relation for the classical decay instability. Hence, 
all conclusions valid for the classical case can be immediately extended to the quantum case. 
In particular, for Q ^> so that the cubic term can be neglected in (|36|). and for uo = ij, 
we obtain the growth rate 

K\E \ 

7 = -^ztT ■ (37) 



"2to 

In all situations, the discriminant 2 - of the cubic equation (}36|) is positive and there are 
one real and two complex conjugate solutions for this equation, one of which is necessarily 
unstable. 

To conclude, there is a formal similarity between the classical and quantum decay insta- 
bilities. The only differences remain in the dispersion relations ()27H29|) . for the quantum 
Langmuir and ion-acoustic modes. The quantum dispersion relations, however, produces a 
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saturation effect not present in the classical case (see Fig. 1). Combining (|29~j) and (J37|) . we 
obtain 



for which a maximum value ^ m ax — \Eq\/\/2H is obtained for K — > oo. This is to be 

compared with the classical case (H = 0) where 7 grows with no bound as K increases. 
Even if the quantum effects do not imply stabilization, they limit the instability to a fixed 
maximum growth rate. 

For dense astrophysical hydrogen plasmas^- where n ~ 10 32 m~ 3 and T e ~ 10 5 K, we 
obtain H ~ 1 and 7 max = |£?o|/v2. For laser hydrogen plasmas, however, to the best of 
our knowledge, such high densities are not yet attainable. For today's typical values^ of 
no ~ 10 28 m~ 3 and T ~ 10^K, we obtain a modest value of H ~ 1CT 2 . The smallness 
of the quantum effects for modulational instabilities in laser plasmas follows from the fact 
that, in this respect, the relevant parameter H is defined as the ratio of the ion plasmon 
energy to the electron thermal energy. If the pertinent quantum parameter were the ratio 
between the electron plasmon energy and the electron thermal energy, as in the quantum 
two-stream instability, 6 the lower mass of the electrons would increase significantly H . For 
laser plasmas with the same typical values as before, we would have H ~ 0.4. Of course, 
we are not saying that quantum effects are irrelevant for laser plasmas: they can show up 
when the Landau length becomes comparable to the de Broglie wavelength, in which case 
huj c ~ KBT e , where uj c is the cyclotron frequency associated with the laser field, or when the 
electrons are degenerated.— 

B. Four-wave instability 

The general instability involving the interaction between one single finite-amplitude 
Langmuir wave, two other Langmuir waves and one ion-acoustic wave can be obtained by 
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(38) 



V2(l + H 2 K 2 y/* ' 




choosing 18 



E(x,t) = E exp(— iu> t + ik x) + E + exp[— i(u + uS)t + i(k + k)x] 





(39) 
(40) 
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where the amplitudes E + , E_ and n are all first order quantities. We choose the equilibrium 
solution E(x,t) = E exp(—iu t + iko%), n(x,t) = to satisfy Eqs. (j23H24j) with E real. 
This implies the relation 

cu = k 2 + H 2 k 4 . (41) 

The last equation can be found also from the dispersion relation for quantum Langmuir 
modes [equation (jHJ)] taking u e = and the recalling of variables. In fact, u e is absorbed in 
the electric field through (|14|). In conclusion, (|41|) shows a quantum Langmuir mode. 
The forms fl39M0|) when inserted in Eq. (J23J) yield 

(lo + u)E + - (k + k) 2 E + = nE + H 2 (k + k) 4 E + , (42) 
(lu - co*)E^ - (k - k) 2 E^ = n*E + H 2 (k - k) 4 E^ . (43) 

When combined, Eqs. (|24p. (j41M2|) and the complex conjugate of ()43|) give the following 
dispersion relation: 

D s D 1 D 2 = k 2 E 2 (D 1 + D 2 ), (44) 

in which 

D s = u 2 — k 2 — H 2 k 4 , (45) 

D 1 = u - k 2 - 2k k - H 2 {k 4 + Ak k 3 + 6k^k 2 + 4JfejjJfe) , (46) 

D 2 = -co - k 2 + 2k k - H 2 (k 4 - Ak k 3 + Qk^k 2 - Ak 3 k) . (47) 

Notice that the limit H — > recovers the classical dispersion relation for the four-wave 
interaction.^ 

The dispersion relation (|44|) is a fourth order polynomial in uj that, in general, can only 
be analyzed numerically. However, the simplest case u> = k = can be investigated 
analytically. For a purely growing instability to = the dispersion relation (|4*4*)l becomes 

[ 7 2 + k 2 + H 2 k 4^ [7 2 + (jfea + R 2 k ^ = 2k *E 2 (k 2 + H 2 k 4 ) . (48) 
Solving for j 2 , there follows two roots, 

7 2 = ~-{k 2 + H 2 k 4 ){l + k 2 + H 2 k 4 ) (49) 
± l -(k 2 + H 2 k 4 ) l/2 [(k 2 + H 2 k 4 )(l - k 2 - H 2 k 4 ) 2 + %k 2 E 2 ] 1 ' 2 , 
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one of them necessarily negative (stable mode). The other root is positive (unstable mode) 
provided 

El > y(l + H 2 k 2 ) 2 . (50) 

This recovers the classical instability condition for the four- wave interaction^ when H — > 0. 
However, there is a new, quantum effect of instability suppression provided 

H>>_^±. (51) 

This is in agreement with the overall stabilization that quantum diffraction effects produce 
in high density plasmas. 6,7 - In fact, for sufficiently large H there is no transfer of energy from 
the original quantum Langmuir mode to the two new quantum Langmuir modes and to the 
quantum ion-acoustic mode. 

Let us consider in more detail the potentially unstable mode described by the positive root 
in flU). In Fig. 2, we show 7 2 versus k 2 for H = 0, H = 0.5 and H = 0.9, with E = 0.5. The 
instability region [ r y 2 > 0) in fc-space becomes narrower for bigger H. Also, the maximum 7 2 
becomes smaller the larger the quantum effects. This latter result is analytically supported 
by an expansion of the positive root of Eq. (J4"9|) . The wavenumber for maximum growth 
rate, k max can be calculated, in a 0(k 5 ) approximation, by expanding Eq. (jlHJ) up to fifth 
order in k. This perturbation analysis is interesting since dj/dk = is not soluble in closed 
form if we use the exact expression for ()49|) . The result of the expansion procedure is shown 
in Fig. 3, where the wavenumber k max for maximum growth rate when E = 0.5 is shown as 
a function of H . (Notice the extended domain of the function, beyond the reasonable limit 
of H = 1.) Using this k max we obtain a somewhat complicated expression which can be used 
to calculate the associated growth rate 7. Using a computer algebra program, we can easily 
obtain ■jmax = Eq(1 — El — 0.87H 2 ), a result valid up to O(Eq). This is an approximate 
equation showing that quantum effects produce stabilization. The approximations adopted 
are justified in view of our assumptions of long wavelengths and weak turbulence (small 
electric field amplitudes). For dense astrophysical plasmas with H ~ 1, as in the decay 
instability case, we would get 7 max = _Eq(0.13 — El), a significant difference in comparison 
with the classical case where ^ m ax = El(l — El). 

In order to further assess the role of quantum effects in the four-wave interaction process, 
we performed a numerical study of f)44j) for general k 7^ 0. Figure 4 displays the real 
(solid lines) and imaginary (dashed lines) parts of uj as a function of k. Both uncoupled 
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(i.e., Eq ~ 0) and coupled cases are considered, for three different values of H. Due to the 
symmetry (k,u) <-> {—k, —uj) =>- D\ <-> D2, of the dispersion relation ([44)1. we consider only 
positive values of the wavenumber, around the overlay region of the branches D s and D 2 , 
where instability occurs. In the uncoupled case, k = 2k is a root of D 2 when uj = 0, for 
both classical and quantum cases. Also, the plots of D s and D 2 branches touch each other at 
isolated points while, when Eq 7^ 0, overlay occurs for a whole finite interval of k, signalizing 
wave instability. The first column of plots shows that, for a fixed ko, both uncoupled curves 
raise with H, implying reduction of the interval in k where instability settles down. This 
can be checked against the corresponding figures in the second and third columns, where a 
contraction of the unstable interval is clearly seen. 

Denote the unstable interval in k by Ik = {k a , For higher pump energy Eq, the third 
column of Fig. 4 shows an overall contraction of Ik- This results from the gradual shift of k a 
to the right and k b to the left, due to the quantum effects. For the relevant range of values 
< H < 1, less severe attenuations occur for the maximum growth rate, compared to those 
found for the unstable interval in k. Thus, the numerical results show that the quantum 
effect inhibits the spreading of energy among different modes. In fact, assume that for a 
specific k, Nj = (k^ — k a )/k represents a first estimation for the number of active modes 
at the beginning of the process. Then, the contraction of Ik implies that the Langmuir 
fluctuations in quantum plasmas might represent more coherent configurations, i.e., having 
less effective modes when compared to the corresponding classical situation, an issue to be 
checked by a direct numerical simulation and, possibly, by an experiment. 



V. NONLINEAR ANALYSIS AND OPEN QUESTIONS 

An important regime of the classical Zakharov equations concerns its static limit. In this 
case, the classical Zakharov system do possess soliton solutions described by a nonlinear 
Schodinger equation— The procedure for the static limit of the quantum Zakharov equations 
considers the approximation d 2 n/dt 2 ~ in ([24)1 . This gives immediately 

Equation inserted in Eq. (J23J), yields 

dE d 2 E o ( E B 2 n\ . . 
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In the classical limit H — > 0, the right-hand side of Eq. (J53*)l vanishes and we recover the 
nonlinear Schrodinger equation with its soliton solutions. In the quantum case, however, 
equations (|52H53|) form a coupled, nonlinear system. We have not been able to find localized, 
analytical solutions for this system. In fact, the usual reduction procedure of searching for 
solutions in the form 



E = F(x — Mt) exp(i[k(x - ut) + 5]) , n = G(x - Mt) 



(54) 



for real F, G, k, M, u and S produces a complicated fourth-order system of coupled, nonlinear 
equations. The existence of soliton solutions for this system remains an open question. It 
seems that a numerical analysis could help in this respect but we believe that this issue 
should be more appropriately treated in a future work. 

Another avenue in nonlinear studies of the quantum Zakharov equation concerns its 
simultaneous semiclassical and static limit. Substituting (J52j) into (|53j) and retaining only 
terms up to 0(H 2 ) produces the decoupled equation 



.dE | d 2 E | 2£ 
dt dx 2 



H 2 



'd 4 E 



E 



d 2 \E\ 



(55) 



v dx 4 dx 2 J 

Equation (|B3j) can be used to study perturbations of the classical NLS soliton solutions. 
The terms proportional to H 2 , in Eq. (|55jl . will probably modify the dispersion-nonlinearity 
equilibrium, which is the ultimate responsible for the soliton existence. 

More formal aspects of the Zakharov equations have to do with its variational formulation 
and the associated Noether currents.— In particular the quantum Zakharov equations pre- 
serve the number of plasmons / \E\ 2 dx of the high frequency electric field, as a consequence 
of the associated conservation law 



dp dJ 
dt dx 







where E(x,t) = A(x, t) exp(i9(x, £)), with A = A(x,t) and 9 
phase functions, and 



P 
J 



A 2 



2 A 



,de 

dx 



2H 2 



A 1 



,d 3 9 

dx 3 



+ 2A 



dA d 2 9 
dx dx 2 



(56) 

(x, t) real amplitude and 

(57) 
(58) 



2 A 



'dfT 
dx 



'dA 
dx 



89 4A d2Ad9 
dx dx 2 dx 
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Notice the extra contribution proportional to H 2 to the plasmons current. The conservation 
law ()56|) comes from the imaginary part of (|23|) and hence contains no contribution from 
n. A proper formulation of the remaining conservation laws (momentum and energy) of 
the system is an open question to be tackled, preferably in accordance with symmetry 
principles of an associated action functional. Other important issues concern the search 
for coherent solutions of the quantum Zakharov equations, namely quantum solitons and 
quantum cavitons. 

Still another issue related to the nonlinear analysis of Eqs. (|23H24|) concerns thermaliza- 
tion and recurrence. For periodic boundary conditions, the classical NLS does not exhibit 
thermalization and, therefore, is generically recurrent . 22 ■ 23 The classical procedure to address 
such questions is based on estimations for the number of active modes Na, from the Rayleigh 
quotient. An upper bound estimation for this number is provided by two invariants: the 
number of plasmons and a momentum-like invariant which, in our case, is not yet known. 
For classical regimes, numerical simulations show that the conclusions can, in general, be 
extended to the non-integrable Zakharov system, when considered as a perturbation of the 
NLS regime.— In fact, it has been shown that, at least for some period of time, the constancy 
of the momentum-like quantity is approximately satisfied. Moreover, numerical simulations 
show that the elementary estimation presented in the last section, i.e., Na ~ Ni, can yield 
quite good results when applied to the full Zakharov equations. Under this viewpoint, the 
contraction of the /c-unstable interval due to H ^ 0, verified in subsection (jIV Bj) . suggests 
that the distribution of energy is less intense in quantum plasmas when compared with the 
classical case. Therefore, quantum effect would favor recurrence in Langmuir modulational 
regimes. 

To finalize, we can derive some exact solutions for the quantum Zakharov equations ([231 
l2"4*j) if we consider pure ion-sound waves obtained by taking E = 0. With zero electric field, 
the density perturbation satisfies the undriven equation 

d 2 n d 2 n - d 4 n „ 

^-a? + ff V = °- (59) 

This linear fourth-order evolution equation was investigated using the method of Lie 
symmetries^ and we found time and space translation symmetries, as well as a scale symme- 
try resulting from the linearity. The H 2 term breaks down the Lorentz invariance endowed 
by the classical model for pure ion-sound waves, so that arbitrary waves travelling at the 
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ion-sound velocity can not be constructed. Nevertheless, exact solutions for Eq. (|59|) can be 
found supposing n = n(x — ct), for constant c and for n a function to be determined. For 
c 2 > 1, corresponding to supersonic flow, and disregarding an integration constant associated 
to non-bound solutions, we get periodic solutions of the form 



where a, b and 5 are numerical constants. This similarity solution is an arbitrary amplitude 
solution. Notice that quantum effects increase the spatial frequency of oscillations in the 
reference frame of the travelling wave. 

VI. CONCLUSION 

We obtained a general model to analyze the coupling between Langmuir waves and ion- 
acoustic waves, in a quantum setting. The model was shown to be appropriate to the 
four-wave interaction and quantum effects have been shown to provide stabilization of a 
classically unstable mode. In the case of the decay instability, a formal similarity with the 
classical case is identified, except for small differences in the dispersion relations, representing 
quantum corrections. We also identified a dimensionless quantum parameter given by the 
ratio of the ion plasmon and electron thermal energies. As pointed out before, this quantum 
parameter may not be small, at least for dense plasmas. 

The consequences of our results on todays laboratory or technological plasmas are not yet 
fully assessed since, for present conditions, H <C 1 in these applications. However quantum 
effects may imply important consequences in the behavior of high density astrophysical 
plasmas, where H ~ 1 is easily found. In this case, as we pointed out, quantum effects cause 
an overall reduction in the wave-wave interaction level. Specifically and in contrast to the 
classical case, the decay instability growth rate is bounded for large wavenumbers. Growth 
rate reduction also occurs for the four-wave interaction. Besides, suppression is also verified 
in the length of the unstable spectral range, implying spectral focusing, i.e., a restriction on 
the range of possible unstable wave-numbers. This focusing effect may extend to quite long 
periods of time, indicating that the recurrence properties verified in the classical Zakharov 
equation are enhanced by the quantum effects. 

A number of open questions remains to be addressed. Of course, a complete analysis of the 




(60) 
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linear dispersion relation of the quantum Zakharov system have to be done. This may require 
a full three-dimensional treatment, with the inclusion of electromagnetic coupling between 
Langmuir and ion-acoustic modes. An additional important point are the nonlinear effects, 
some of them briefly discussed in Section V, which may deserve a more careful scrutiny. To 
conclude, the huge amount of physical and mathematical aspects already assessed in the 
classical Zakharov equations certainly have quantum counterparts which ask for an equally 
careful investigation. 
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FIGURE CAPTIONS 



FIG. 1. Growth rate of the decay instability for Eq = 0.5 and < H < 1.5, as indicated. 
Notice the quick saturation effect for H > 0. 



FIG. 2. 7 2 as a function of k 2 for the positive root in the dispersion relation ()49|) for 
the four-wave instability. We have Eq = 0.5, H = (full line), H = 0.5 (dashed line) and 
H = 0.9 (dotted line). 



FIG. 3. Wave-number k max for maximum growth rate of the four-wave instability, as a 
function of H, calculated to 0(k 5 ) and E = 0.5. 



FIG. 4. Real (solid lines) and imaginary (dashed lines) components of the frequency uj 
as a function of k for uncoupled (frames al, a2 and a3) and coupled cases (frames bl to c3). 
From top to bottom, H = 0, H = 0.5 and H = 0.9, respectively. From left to right, E = 0, 
Eq = 0.5 and Eq = 0.5. For the first and second columns fco — 0.5; for the third column, 
k = 0.75. In the first frame, Di (i = s, 1, 2) indicate the various branches of (}4*4"|) . A similar 
labelling applies to all the frames. 
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